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INTRODUCTION 

Use  of  gadolinium-based  magnetic  resonance  (MR)  contrast  agents  (gadopentetate 
dimeglumine)  is  known  to  offer  exceptionally  high  sensitivity  for  the  detection  of  breast 
cancers  [1-6].  The  specificity  of  gadolinium-enhanced  MR  mammography  to  distinguish 
benign  from  malignant  tissues,  however,  has  been  limited  by  the  fact  that  some  normal 
and  benign  masses  also  exhibit  substantial  contrast  enhancement  [3,5,6].  Alternatively, 
the  "temporal  rate"  of  gadolinium  uptake  is  being  studied  as  a  potential  complimentary 
discriminator  of  benign  from  malignant  breast  neoplasm  [7-10].  Malignancies  typically 
exhibit  relatively  rapid  contrast  enhancement  by  virtue  of  their  anomalous  vascularity  in 
terms  of  high  capillary  density  and/or  highly  permeable  capillaries.  Non-invasive 
methods  to  detect  and/or  discriminate  between  normal,  benign,  pre-malignant  and 
malignant  tissues  based  on  vascular  alterations  would  have  significant  diagnostic  value, 
particularly  in  cases  that  are  a  challenge  for  conventional  mamography  such  as  in  dense 
breasts. 

Acquisition  and  processing  methods  to  quantify  MRI  contrast  agent  enhancement  rates, 
or  temporal  properties  in  general,  have  been  widely  variable  across  studies  [7-10].  For 
accurate  characterization  of  contrast  uptake  features,  the  most  critical  parameters  include: 
temporal  resolution,  spatial  resolution,  volume  of  tissue  coverage,  and  the  method  to 
quantify  signal  changes.  Information  regarding  “rate  of  enhancement”  is  obscured  if  the 
temporal  sampling  is  inadequate  to  record  details  of  the  lesion  enhancement  curve  such  as 
the  point  contrast  onset  and  maximal  slope.  Problems  also  arise  if  the  imaged  volume 
only  partially  includes  an  enhancing  lesion  or  does  not  include  all  lesions  in  the  breast. 
These  issues  arise,  in  part,  due  to  the  well  known  trade  off  between  volume  coverage, 
spatial  resolution,  and  temporal  resolution. 

The  central  objectives  of  this  project  are  technical  refinement  and  clinical 
implementation  of  new  dynamic  contrast-enhanced  MRI  methods  designed  to  provide 
full  volume  coverage  of  both  breasts  at  good  spatial  and  temporal  resolution  for  detailed 
quantification  of  dynamic  contrast  enhancement  properties  [10].  This  approach  involves 
specialized  software  for  data  acquisition  (i.e.  MRI  pulse  sequence),  image  reconstruction, 
and  enhancement  curve  data  reduction  of  gadolinium-enhanced  tissues.  The 
methodology  allows  dynamic  scans  through  32  sections  of  both  breasts  acquired  at  10- 
13second  temporal  resolution  -  this  is  4-fold  faster  than  the  conunercially  available 
sequence.  Additional  post-processing,  referred  to  as  “viewsharing”,  further  reduces 
temporal  sampling  to  5-6seconds.  This  technique  is  being  applied  to  study  100  women  in 
whom  a  breast  abnormality  detected  by  mammography,  physical  exam,  or  ultrasound  is 
to  be  biopsied. 
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METHODS 

Acquisition  Method 

Breast  MR  exams  are  performed  on  l.STeslta  General  Electric  Signa  systems  using  a 
conunercial  (GE)  bilateral  phased-array  breast  coil  with  insert  pads  for  mild 
anterio/posterior  (A/P)  compression.  After  a  quick  locator  scan,  conventional  axial  "T2- 
weighted"  fast-spin-echo  is  performed.  This  sequence  is  primarily  used  to  identify  cysts. 
Axial  pre-  and  post-gadolinium  "T1 -weighted"  spin-echo  bracket  the  "dynamic  scan" 
described  below. 

A  volumetric  rf-spoiled  gradient  recalled  echo  (3D  SPGR)  sequence  is  used  for  the 
dynamic  acquisition  since  it  (a)  offers  excellent  sensitivity  to  gadolinium-induced 
changes,  (b)  produces  high  signal-to-noise  images,  and  (c)  provides  full  volume  coverage 
of  both  breasts  with  contiguous  image  sections.  Default  parameter  settings  include: 
Coronal  planes;  TR=10ms;  TE=4.5ms;  flip=40O;  FOV=30cm  right-left  (R/L),  15cm 
superior-inferior  (S/I),  3mm  slice  thickness;  acquisition  matrix=256  (R/L),  128  (S/I),  32 
slices  (A/P). 

Scan  time  and  spatial  resolution  are  proportional  to  the  acquisition  matrix  size. 

Moreover,  dynamic  changes  (due  to  gadolinium)  are  most  heavily  concentrated  near  the 
center  of  the  3D  Fourier-encoded  dataset  commonly  referred  to  as  “k-space”.  Therefore, 
to  successfully  record  most  of  the  dynamic  contrast  information  and  greatly  reduce  scan 
times  for  high  temporal  resolution  (i.e.  image  set  every  10  seconds),  only  the  central  k- 
space  data  are  acquired  during  the  gadolinium  injection.  The  lost  spatial  resolution  (since 
peripheral  k-space  data  were  not  acquired)  is  partially  restored  by  reuse  of  previously 
acquired  (non-dynamic)  peripheral  k-space  data.  This  process  of  splicing  high  temporal 
resolution  central  k-space  data  blocks  with  high  spatial  resolution  peripheral  k-space  data 
is  referred  to  as  the  “keyhole”  method,  and  has  been  extended  to  3D  for  this  study  [10]. 

An  additional  post-processing  technique,  referred  to  as  “viewsharing”,  is  applied  to 
effectively  double  the  temporal  sampling  rate  (i.e.  image  set  every  5  seconds). 
Viewsharing  is  the  synthesis  of  data  blocks  at  each  intermediate  time  point  by 
combination  of  the  latter  half  of  a  given  data  keyhole  block  with  the  initial  half  of  the 
subsequent  data  block.  As  with  each  original  keyhole  data  block,  synthetic  viewshared 
blocks  are  submitted  for  keyhole  reconstruction.  In  this  study,  a  series  of  21  keyhole 
blocks  are  acquired  over  a  4-5minute  period.  The  standard  dose  of  gadolinium  contrast 
material  (O.lmM/kg)  is  administered  as  a  bolus  30  seconds  into  the  dynamic  scan. 
Viewsharing  interleaves  additional  21  points  for  a  42-point  total  time  series. 

Processing  Methods 

Raw  data  are  networked  to  a  SUN  SPARC  10  for  offline  3D  keyhole/viewshare 
processing.  This  involves  some  corrections  and  insertion  of  each  keyhole  dataset  into  the 
full  matrix  reference  dataset,  then  3D  Fourier  transform  to  yield  32  slices  at  each 
timepoint.  The  last  pass  dataset  is  used  to  correct  an  enhancement  amplitude  error  that 
results  from  partial  k-space  acquisition  [10].  In  addition,  a  motion  correction  algorithm  is 
applied  to  independently  register  each  breast  at  each  timepoint  to  the  reference  data  set 
(i.e.  pre-injection).  The  correction  algorithm  assumes  the  object  is  limited  to  simple 
^splacment  in  3-dimensions.  Motion-corrected  keyhole  reconstructed  images  and 
"reference-subtracted"  images  are  viewed  cinegraphically  on  the  workstation.  Rapidly 
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enhancing  lesions  are  most  conspicous  when  viewed  in  this  format,  thus  assisting  the  user 
in  region  of  interest  (ROI)  definition. 

Lesion  ROI  enhancement  time-course  curves  are  automatically  fit  to  mono-  and  double¬ 
exponential  saturation,  and  saturation-decay  exponetial  models.  These  mathematical 
functions  are  derived  from  compartmental  models  that  describe  the  temporal  distribution 
of  contrast  media  between  tissue  vascular  and  extracellular  spaces.  These  models,  albeit 
simplistic,  suggest  that  capillary  permeability  and  surface  area  (P*S  product)  are  most 
relavant  to  early  enhancement  rate.  Enhancement  rate  is  defined  as  [fitted  enhancement 
time,  tau]*\  Also,  overall  enhancement  amplitude  (i.e.  degree  of  signal  change)  reflects 
size  of  extracellular  space.  Fitted  enhacement  amplitude  and  temporal  properties  are 
automatically  recorded  on  disk  for  each  ROI.  Image  excerps,  contrast  enhancement 
curve,  and  fitted  results  for  an  invasive  carcinoma  are  illustrated  in  Figure  1. 


Invasive  Ductal  Ca 


Figure  1.  Excerpts  from  an  enhancement  time  series  of  a  invasive  ductal  carcinoma.  This 
lesion  exhibits  a  high  enhancement  rate  and  "washout"  (signal  decay);  both  signs  are 
strongly  suggestive  of  cancer. 
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Statement  of  Work  Tasks  which  relate  to  technical  developments  are: 

Task  2:  F.nhancement  rate-constant  determination  -  each  patient 

Each  patient  will  be  studied  by  a  new  dynamic  MRI  sequence  during  bolus 
administration  of  gadolinium-DTPA  contrast.  The  data  acquisition  and  processing 
sequence  will  generate  a  temporal  series  of  32  data  sets,  where  each  set  represents  32 
tomographic  slices  through  both  breasts,  acquired  at  a  rate  of  1  set  per  10-13  seconds. 
Existing  specialized  interactive  software  will  be  employed  to  retrospectively  view 
each  tomographic  slice  cinegraphically  to  identify  gadolinium  enhanced  tissues. 
Enhancement  rate-constants  will  be  quantified  by  numerical  fit  to  signal  change 
versus  time  using  a  exponential  saturation  model. 

Task  3:  Software  development:  View  Sharing  &  Automated  Fitting  -  1st  year 

View  sharing  is  a  post-processing  concept  designed  to  effectively  double  the 
dynamic  sampling  rate  using  existing  data.  Software  to  execute  and  test  this  concept 
will  be  written  for  simulated  and  actual  MRI  data.  Pending  satisfactory  performance 
of  view  sharing,  archived  data  sets  will  be  reprocessed,  and  subsequent  data  sets  will 
be  processed  with  view  sharing  thus  yielding  64  data  sets  at  6sec  dynamic  temporal 
resolution.  Addition  development  work  will  be  directed  toward  export  of  datasets 
from  workstation  back  into  MRI  system  database  for  more  portable  viewing. 

These  tasks  have  been  completed  with  the  exception  of  exporting  datasets  back  into  the 
MRI  database  for  more  portable  viewing.  Archived  datasets  have  been  processed  with 
and  without  viewsharing,  as  well  as,  with  and  without  motion  correction.  Thesis  work  of 
Sumati  Krishnan,  the  biomedical  engineering  student  supported  by  this  grant,  involves 
comparison  of  these  methods.  We  have  elected  to  retain  dynamic  breast  image  datasets 
on  the  workstation  since  transferal  and  storage  of  large  image  sets  (>2000  images  per 
case)  is  cumbersome.  Instead,  we  have  added  high-quality  hardcopy  capability  and  a 
2.6Gbyte  magneto-optical  drive  to  the  SPARC  10  workstation. 

Clinical  Method 

Recruited  patients  are  referred  from  the  Breast  Imaging  Section  of  the  Breast  Care  Center 
at  the  University  of  Michigan  Hospital.  Entrance  criteria  includes  any  patient  with  a 
breast  abnormality  detected  by  mammography,  physical  examination,  or  ultrasound 
which  will  be  biopsied.  Prior  to  MRI  study,  all  patients  have  a  routine  clinical  evaluation 
which  includes,  as  a  minimum,  mammography  and  breast  physical  examination. 

Three  radiologists  independently  review  the  data  from  the  following  "modalities"  in  the 
order  of:  (a)  conventional  mammography  images  (Xray,  ultrasound)  and  reports;  (b) 
conventional  MRI  including  pre-  and  post-gadolinium;  and  (c)  subtracted  MR  images 
(pre-gad  subtracted  from  post-gad).  Readers  grade  lesions  on  five  point  scales  of 
suspicion  (benign  to  malignant)  and  visibility  (obvious  to  subtle). 


Statement  of  Work  Tasks  which  relate  to  clinical  methods  are: 

Task  1:  Patient  recruitment  -  uniform  over  the  four  years  of  the  study 
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One-hundred  breast  patients  will  be  studied  during  this  four-year  prospective 
imaging  trial.  Volunteer  patients  will  be  recruited  from  the  Breast  Imaging  Section  of 
the  Breast  Care  Center  at  the  University  of  Michigan  Hospital.  Entrance  criteria  will 
include  patients  with  a  breast  abnormality  detected  by  mammography,  physical  exam, 
or  ultrasound  which  will  be  biopsied.  Malignancy,  grade,  type  and  lymph  node 
involvement  will  be  determined  for  each  patient  by  pathologic  reading. 
Approximately  30  malignant  and  70  benign  lesions  are  expected  from  this  population. 

Task  4:  Statistical  summary  -  final  two  months  of  study 

Quality  of  curve  fitting  and  appropriateness  of  selected  model  will  be  statistically 
assessed  on  individual  patient  data.  This  will  be  performed  throughout  the  study,  as 
will  assembly  of  other  diagnoses  and  pathological  outcomes  into  a  database.  Final 
inter-patient  statistical  summaries  will  occur  in  the  last  few  months  of  the  study. 
Quantitative  enhancement  rate  constants,  conventional  pre-biopsy  diagnosis,  and 
pathologic  characterization  will  be  submitted  for  statistical  testing  of  stated 
hypotheses  regarding  diagnostic  sensitivity,  specificity,  positive  and  negative 
predictive  values  of  enhancement  rate  constants  relative  to  conventional  pre-biopsy 
diagnosis.  Classification  power  for  tissue  group  type,  and  correlation  with  tumor 
grade,  size,  invasiveness,  and  lymph  node  involvement  will  be  determined. 

As  of  Oct  '97, 65  patients  have  been  recruited  into  the  study.  While  this  is  lower  than  the 
projected  number,  we  do  not  anticipate  difficulty  meeting  our  target  of  100  patients. 

Task  4  will  be  formally  addressed  towards  the  end  of  the  study;  although  the  UM 
Comprehensive  Cancer  Center,  Biostatics  Core  Director,  KyungMann  Kim,  Ph.D.,  has 
identified  graduate  student  support  for  this  project. 

RESULTS 

Spatial  correlation  between  pathology,  radiologic  follow-up  (>2years)  and  ROIs 
quantified  by  MRI  was  achieved  for  61  lesions/tissues  as  summarized  in  the  Table  below. 

Table 


Mean  ±  StdDev  of  MRI-derived  Contrast  Enhancement  Amplitude  and  Rate  Properties 


Enhancement 

Enhancement 

Enhancement 

Tissue  Type 

Rate  (~  P*S) 

[10'^  sec‘] 

Amplitude 

Amplitude  x 

[%  Signal  Change] 

Enhancement  Rate 

Benign  (N=38) 

16  ±10 

146  ±99 

3±3 

Malignant  (N=23) 

58  ±46 

174  ±90 

10  ±10 

p  value 

<  0.001 

N.S. 

<0.01 

Scatterplots  for  individual  tissues  of:  (a)  enhancement  rate,  which  is  primarily  dependent 
on  PS  product  and  (b)  enhancement  amplitude  are  provided  in  Figure  2. 


Page  6 


Grant  Number:  DAMD17-94-J-4381 


ff 


(a)  Enhancement  rate 


^Fibrocy  ch 
H  Fibroadenoma 
A  Benign 
X  Nippb 

■  LCB 
•  DCB 

<^iivDuctalCa 
APapilkry  Ca 
#MeduIlary  Ca 

■  iiv  Lobular  Ca 


ec 

fl 

ed 


ed 

a 

.2f 

0) 

TS 

0 

ft 

S 

<< 

• 

ja 

B 

u 


450 

400 

350 

300 

250 

200 

150 

100 

50 

0 


♦ 


^Fibrocy  ch 

■  Fibroadenoma 
^Benign 

X  Nipple 

eLCB 

®DCB 

♦  hvDuctalCa 
^Papillary  Ca 

•  Medullary  Ca 

■  Iiv  Lobular  Ca 

10  11 


Figure  2.  Scatterplot  summary  of  quantitative  enhancement  properties:  (a)  enhancement 
rate,  and  (b)  enhancement  amplitude  sorted  by  tissue  types. 
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CONCLUSIONS 

Technical  objectives  of  high  temporal  sampling  and  good  spatial  resolution  for 
quantitative  description  of  breast  lesion  enhancement  properties  are  well  satisfied  by  this 
methodology.  Motion  correction  prior  to  keyhole  reconstruction  provides  a  qualitative 
improvement  by  reducing  motion  artifact  in  subtracted  images.  Analysis  of  enhancement 
properties  indicates  the  rate  of  enhancement  is  a  better  discriminator  of  malignancy  than 
the  amount  (amplitude)  of  enhancement.  Invasive  lobular  carcinoma  (N=3)  and  a 
minority  of  invasive  ductal  carcinomas,  however,  did  not  exhibit  the  typical  high  rate  of 
enhancement  of  other  cancers  resulting  in  an  overlap  between  benign  and  malignant 
tissues.  Nevertheless,  high  sensitivity  of  dynamic  contrast-enhanced  MRI  to  lesions,  and 
the  ability  to  detect  additional  lesions,  even  in  dense  breasts,  suggests  these  methods  may 
have  value  in  the  following  clinical  situations:  (1)  Survey  breasts  for  additional  foci  prior 
to  surgical  biopsy  of  lesions  detected  by  conventional  means,  particularly  in  dense  tissues 
or  those  obscured  by  implants.  (2)  Means  of  increased  surveillance  of  women  who  "test 
positive"  by  future  serologic  or  genetic  breast  cancer  screens. 
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Djmamic  Three-dimensional  Imaging 
with  Partial  K-Space  Sampling:  Initial 
Application  for  Gadolinium-enhanced  Rate 
Characterization  of  Breast  Lesions^ 


PURPOSE:  To  evaluate  a  method  to 
monitor  gadolinium  enhancement  pat¬ 
terns  at  magnetic  resonance  (MR)  imag¬ 
ing  with  high  temporal  resolution  and 
full  coverage  through  both  breasts. 

MATERIALS  AND  METHODS:  In  12 
patients  with  13  masses,  including 
nine  carcinoma,  nonenhanced  three- 
dimensional  MR  imaging  was  per¬ 
formed  with  full-matrix  resolution. 

At  dynamic  imaging,  32  serial  passes 
were  made  during  bolus  administra¬ 
tion  of  contrast  material,  and  tempo¬ 
ral  resolution  was  reduced  to  12  sec¬ 
onds  by  collecting  the  central  (low 
spatial  frequency)  32  x  16  or  16  x  16 
phase-encode  views.  Full-matrix  dy¬ 
namic  images  were  reconstructed  by 
complementing  central  phase-encode 
data  with  precontrast  data  from  pe¬ 
ripheral  high-spatial-frequency  views. 

RESULTS:  Results  at  time-course 
analysis  with  a  mono-exponential 
saturation  model  indicated  malignant 
lesions  tend  to  show  rapid  (<60  sec¬ 
onds)  contrast  change  relative  to  be¬ 
nign  masses  and  normal  tissues.  One 
cancer  displayed  an  exceptionally 
slow  contrast  change  (260  seconds). 

CONCLUSION:  The  technical  objec¬ 
tives  of  full  tissue  coverage,  rapid 
temporal  sampling,  and  quantifica¬ 
tion  of  enhancement  curves  are  met 
with  this  method  for  certain  lesions 
(>5  mm  in  largest  diameter). 

Index  terms:  Breast,  MR,  00.12143  •  Breast 
neoplasms,  MR,  00.12143  •  Magnetic  resonance 
(MR),  k-space 
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ARLY  magnetic  resonance  (MR) 
imaging  demonstrated  only  lim¬ 
ited  success  in  assisting  diagnosis  of 
breast  disease  (1-4).  Initially,  only  tis¬ 
sue-inherent  differences  in  nuclear 
MR  relaxation  times,  T1  and  T2,  were 
the  primary  sources  of  contrast.  Al¬ 
though  these  contrasts  are  relatively 
high,  it  was  observed  that  they  have 
little  value  in  helping  detection  and 
diagnosis  of  breast  cancer  (1,3,5,6). 
Variable  admixtures  of  intense  fat  sig¬ 
nal  and  inadequate  fat  suppression 
have  also  seriously  hindered  image 
interpretation  and  detection  of  breast 
carcinoma. 

More  recently,  gadolinium-en¬ 
hanced  MR  mammography  has  gen¬ 
erated  considerable  interest  (7-10).  It 
is  generally  accepted  that  gadolinium 
enhanced  MR  imaging  offers  high 
sensitivity  for  the  detection  of  breast 
cancers.  The  specificity  of  gadolinium- 
enhanced  MR  mammography  to  dis¬ 
tinguish  benign  from  malignant  tis¬ 
sues,  however,  has  been  limited  by 
the  fact  that  benign  tissues  can  also 
exhibit  contrast  enhancement  (9,11,12). 

The  rate  of  gadolinium  enhance¬ 
ment  has  also  been  studied  as  a  po¬ 
tential  discriminator  of  benign  and 
malignant  breast  neoplasms  (13-15). 

In  their  initial  summary  of  findings  in 
25  patients  who  underwent  dynamic 
gadolinium-enhanced  MR  imaging, 
Kaiser  and  Zeitler  reported  that  all 
malignancies  displayed  rapid  en¬ 
hancement  (ie,  ~  100%  increase  in 
signal  intensity  within  the  first  2  min¬ 
utes  after  injection  of  0.1  mmol/kg 
gadopentetate  dimeglumine),  whereas 
benign  lesions  demonstrated  less  and 
slower  signal  intensity  change  (13), 


These  authors  have  since  updated 
their  dynamic  gadolinium-enhanced 
MR  imaging  study,  in  which  sensitiv¬ 
ity  and  specificity  for  cancer  detection 
was  greater  than  95%  (16).  In  other 
dynamic  gadolinium-enhanced  MR 
imaging  studies,  higher  false-positive 
rates  were  observed;  thus,  there  is  not 
a  consensus  among  investigators  on 
methods  or  diagnostic  value  of  this 
approach  (11,12,17).  Findings  in  these 
and  other  studies  illustrate  that  fibro¬ 
adenomas  and  proliferative  fibrocystic 
change,  in  particular,  are  sources  of 
false-positive  findings.  Breast  malig¬ 
nancies  typically  exhibit  rapid  en¬ 
hancement  while  benign  lesions  usu¬ 
ally,  although  not  necessarily,  exhibit 
slower  contrast  change.  Malignant 
lesions  may  enhance  slowly,  thus  re¬ 
ducing  sensitivity  (12,18). 

Use  of  various  acquisition  and  en¬ 
hancement  rate  quantification  schemes 
may  contribute  to  discrepancies  in 
findings  (13-15).  In  assessment  of 
enhancement  rate,  important  param¬ 
eters  include  temporal  resolution, 
volume  of  tissue  covered,  and  the 
method  used  to  quantify  signal  inten¬ 
sity  change.  In  particular,  important 
information  on  the  rate  of  enhance¬ 
ment  may  be  obscured  if  the  temporal 
resolution  is  inadequate  to  sample 
rapid  contrast  changes.  In  addition, 
problems  will  obviously  arise  if  the 
dynamic  imaging  volume  includes 
an  enhancing  lesion  only  partially 
or  does  not  include  all  lesions  in  the 
breast.  These  issues  arise  in  part  be¬ 
cause  there  is  an  inevitable  trade-off 
in  standard  MR  imaging  between  im¬ 
aged  volume  and  temporal  resolution. 

MR  imaging  with  methods  such  as 
high  spatial  resolution,  three-dimen¬ 
sional  (3D)  volumes,  suppression  of 
fat  signal,  and  magnetization  transfer 
contrast  have  been  studied  in  an  at- 


Abbreviations:  ROI  =  region  of  interest,  3D  = 
three-dimensional. 
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tempt  to  produce  high-quality  ana¬ 
tomic  images  (9,19).  These  3D  acquisi¬ 
tions  typically  require  several  minutes; 
thus,  dynamic  properties  are  not  readily 
observed.  Conversely,  to  improve 
temporal  sampling  of  enhancement 
curves,  others  have  performed  two- 
dimensional  multisection  MR  imaging 
with  1 -minute  temporal  resolution 
(13)  or  of  one  judiciously  selected  sec¬ 
tion  with  finer  temporal  (15-second) 
resolution  (14).  Unfortunately,  selection 
of  prospective  section(s)  for  dynamic 
MR  imaging  by  means  of  correlation 
with  findings  at  x-ray  mammography, 
ultrasonography  (US),  or  physical  ex¬ 
amination  can  be  extremely  difficult 
in  the  amorphous  anatomy  of  a  vari¬ 
ably  compressed  breast,  except  in 
cases  of  large  tumors.  Moreover,  the 
nonspecific  features  of  some  lesions 
on  precontrast  MR  images  may  fur¬ 
ther  hinder  prospective  selection  of 
appropriate  section(s)  for  dynamic 
imaging. 

The  central  objective  of  this  study 
was  to  obviate  the  trade-off  between 
tissue  coverage  and  temporal  resolu¬ 
tion  seen  in  other  studies  of  dynamic 
gadolinium-enhanced  MR  imaging  of 
the  breast.  The  specific  technical  ob¬ 
jectives  were  development  and  vali¬ 
dation  of  a  pulse  sequence  that  would 
{a)  monitor  gadolinium  enhancement 
on  a  short  time  scale  ( 12  seconds), 

(b)  image  the  full  volume  of  both 
breasts,  (c)  produce  anatomic  images 
with  sufficient  spatial  resolution  and 
contrast  to  allow  correlation  of  find¬ 
ings  with  those  on  images  obtained 
with  other  modalities,  and  (d)  quan¬ 
tify  parameters  related  to  enhance¬ 
ment  rate. 

MATERIALS  AND  METHODS 
Patients 

Twelve  female  patients,  aged  30-62 
years  (mean,  48  years),  were  recruited  for 
this  study.  Eleven  patients  had  one  or 
multiple  masses  identified  at  mammogra¬ 
phy  and/or  physical  examination.  Patho¬ 
logic  confirmation  was  obtained  in  11  of  13 
masses,  including  nine  invasive  carcinoma 
tumors  (largest  diameter  mean,  22  mm; 
range,  12-50  mm)  in  seven  patients  and  in 
the  one  benign  mass  in  each  of  two  pa¬ 
tients  (fibrosis  with  largest  diameter  of  30 
mm  and  a  palpable  benign  mass  with  larg¬ 
est  diameter  of  10  mm).  Of  the  remaining 
two  of  13  masses,  one  was  a  palpable  mass 
with  a  largest  diameter  of  20  mm  that  re¬ 
solved  with  time  in  one  patient  and  was 
therefore  presumed  to  be  benign  and  the 
other  was  a  mass  with  a  largest  diameter 
of  17  mm  that  was  radiographically  consis¬ 
tent  with  a  fibroadenoma  and  was  stable 
for  7  years  in  another  patient.  Cytologic 
findings  were  also  consistent  with  a  fibro¬ 


adenoma.  One  patient  had  no  detectable 
mass  and  was  presumed  to  have  no  breast 
disease.  The  MR  imaging  methods  were 
approved  by  an  institutional  review 
board,  and  informed  written  consent  was 
obtained  from  all  participants. 


Data  Acquisition  for  3D  Keyhole 
MR  Imaging 

Patients  underwent  imaging  on  a  Signa 
imager  (GE  Medical  Systems,  Milwaukee, 
Wis).  Plane  orientation  and  field  of  view 
for  the  dynamic  series  were  chosen  to  al¬ 
low  efficient  imaging  of  the  breasts  with 
high  spatial  resolution  and  minimal  image 
wraparound  and  involvement  of  the 
heart.  Coronal  imaging  was  performed 
through  the  suspended  breasts  with  use 
of  a  rectangular  field  of  view  with  matrix 
dimensions  of  256  (right  to  left)  x  128  (su¬ 
perior  to  inferior)  pixels  as  the  preferred 
geometry.  In  eight  patients  imaging  was 
performed  with  a  dedicated  breast  coil. 

In  four  other  patients  imaging  was  per¬ 
formed  with  the  patient  lying  prone  on 
a  surface  coil,  as  a  result  of  body  habitus, 
and  axial  dynamic  imaging  was  necessary 
in  two  of  these  patients.  The  remainder  of 
the  dynamic  studies  were  performed  with 
coronal  views. 

The  patients  initially  underwent  local¬ 
ization  with  conventional  nonenhanced 
Tl-  and  T2-weighted,  spin-echo  se¬ 
quences.  Conventional  enhanced  spin- 
echo,  Tl-weighted  imaging  was  also  per¬ 
formed.  A  reduced  k  space  or  "'keyhole" 
technique  (20-22),  described  in  detail  be¬ 
low,  was  used  for  dynamic  imaging.  Pre¬ 
contrast  acquisitions,  referred  to  as  refer¬ 
ence  data  sets,  were  performed  with  a  3D 
matrix  of  256  x  128  x  32  sections  and  one 
signal  acquired. 

Reference  imaging  was  followed  imme¬ 
diately  with  dynamic  3D  keyhole  imaging, 
without  changing  the  hardware  settings, 
and  only  the  central  phase-encode  views 
(in  both  phase-encode  directions)  were 
serially  imaged  in  32  passes  during  admin¬ 
istration  of  contrast  material.  In  nine  of 
the  12  patients  studied,  3D  keyhole  dimen¬ 
sions  were  256  x  32  x  16  (ie,  one-eighth 
of  k  space  acquired),  and  in  the  remain¬ 
ing  three  patients  3D  keyhole  dimensions 
were  256  x  16  x  16  (ie,  one-sixteenth  of 
k  space  acquired).  Timing  parameters  in¬ 
cluded  repetition  time  of  24  msec  and 
echo  time  of  4.5  msec  (24/4.5)  for  the 
256  X  32  X  16  keyhole  acquisitions  and 
45/7  for  the  256  x  16  x  16  keyhole  acqui¬ 
sitions.  Temporal  resolution,  defined  as 
the  time  to  complete  each  keyhole  pass, 
was  comparable  in  all  acquisitions  (ie,  12.3 
and  11.5  seconds,  respectively).  Reference 
and  dynamic  keyhole  imaging  were  both 
performed  with  use  of  the  same  radio-fre¬ 
quency-spoiled  sequence  of  gradient-re- 
called  acquisition  in  the  steady  state, 
which  included  a  frequency-selective,  fat- 
saturation  pulse.  Other  imaging  param¬ 
eters  included  a  40°  flip  angle,  one  signal 
acquired,  32-bit  digitization,  and  ±32-kHz 
bandwidth.  The  nominal  right-to-left  field 
of  view  of  300  mm  and  section  thickness 


of  3-6  mm  (mean,  4.5  mm)  were  adjusted 
for  patient  size.  The  position  and  thick¬ 
ness  of  3D  slabs  (equal  to  28  sections)  were 
graphically  prescribed  on  orthogonal  im¬ 
ages  to  encompass  the  volume  of  the 
breasts.  Imaging  times  for  reference  and 
serial  3D  keyhole  acquisitions,  respec¬ 
tively,  were  1  minute  39  seconds  and  6 
minutes  33  seconds  with  repetition  time 
of  24  msec,  and  3  minutes  5  seconds  and  6 
minutes  9  seconds  with  repetition  time  of 
4  msec. 

In  all  cases,  0,1  mmol/kg  of  gadopen- 
tetate  dimeglumine  (Magnevist;  Berlex 
Laboratories,  Wayne,  NJ)  was  manually 
injected  as  a  bolus  over  10-15  seconds  in 
an  antecubital  vein.  The  injection  was 
started  5  minutes  30  seconds  before  the 
end  of  dynamic  imaging.  After  dynamic 
imaging,  the  raw  data  were  permanently 
stored  (requiring  1-2  minutes).  Then  a 
full-matrix  series  was  obtained  with  the 
same  imaging  parameters  used  for  the  ref¬ 
erence  data  set  to  enable  comparison  of 
pre-  and  postcontrast  images. 


Data  Processing 

Time-domain  data  from  precontrast  and 
dynamic  acquisitions  were  transferred  to  a 
computer  (SparclO;  Sun  Microsystems, 
Mountain  View,  Calif)  for  off-line  process¬ 
ing.  Custom  software  tools  for  reconstruc¬ 
tion,  display,  and  region-of-interest  (ROI) 
analysis  were  written  for  processing  (Ad¬ 
vanced  Visual  Systems,  Waltham,  Mass). 
Modules  were  also  written  to  determine 
phase  differences  between  the  early  dy¬ 
namic  passes  (ie,  prior  to  changes  in  con¬ 
trast)  and  the  corresponding  central  time- 
domain  data  of  the  full-matrix  reference 
data  set.  Ideally,  these  should  have  identi¬ 
cal  phase.  Slight  patient  repositioning, 
however,  could  result  in  phase  error  and 
artifact  on  reconstructed  keyhole  images. 
Therefore,  phase  differences  were  deter¬ 
mined  and,  when  necessary,  a  single 
phase-correction  setting  was  applied  to 
all  32  temporal  passes  of  each  study. 

After  phase  correction,  each  pass  was 
substituted  for  the  central-time-domain 
core  of  the  full-matrix  reference  data  set 
and  was  reconstructed  into  a  32-section 
image  set  by  means  of  conventional  3D 
Fourier  transform.  This  process  was  re¬ 
peated  for  all  passes  to  produce  1,024 
anatomic  images  (32  sections  x  32  time 
points).  Subtraction  images  were  also  pro¬ 
duced  by  using  an  early  time  point  as  a 
mask.  Anatomic  and  subtraction  images 
versus  time  were  viewed  side  by  side  cine- 
graphically.  Other  display  formats  de¬ 
signed  for  quick  survey  of  data  sets  in¬ 
cluded  an  image  matrix  of  all  sections 
(anatomic  and  subtracted  images)  at  the 
last  time  point  and  cinegraphic  display 
of  maximum  intensity  projections  of  sub¬ 
tracted  images.  The  maximum-intensity- 
projection  images  were  obtained  through 
the  full  tissue  volume  along  orthogonal 
views  and  were  particularly  beneficial  for 
viewing  changes  in  spatial  and  temporal 
contrast  of  the  lesion(s)  relative  to  the 
heart  and  large  vessels. 
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The  time  course  of  signal  change  in  an 
ROI  of  enhancing  tissues  was  fit  to  a 
mono-exponential  saturation  model  with 
use  of  a  Levenberg-Marquardt  least- 
squares  optimization  algorithm  (23).  Fitted 
parameters  include  contrast  onset,  am¬ 
plitude  of  signal  change,  A;  and  enhance¬ 
ment  time  constant,  t,  used  within  the  fol¬ 
lowing  model: 


AS(t)  = 


A[1  -  for  t  >  to 

I  0,  for  t  <  to. 


(1) 


RESULTS 

Regarding  mass  conspicuity,  con¬ 
ventional  precontrast  spin-echo  Tl- 
and  T2-weighted  images  frequently 
did  not  delineate  location  or  spatial 
extent  of  the  masses.  Visual  compari¬ 
son  of  pre-  versus  postcontrast  images 
(Tl-weighted  spin-echo  and  gradient- 
recalled-echo  images)  was  adequate 
to  identify  nearly  all  (12  of  13)  pal¬ 
pable  and/or  mammographically  de¬ 
tected  masses.  Mass  detection  and 
localization  were  facilitated  by  inspec¬ 
tion  of  anatomic  and  subtracted  cine- 
graphic  images  generated  by  means 
of  the  keyhole  scheme.  Subtraction 
removed  many  image  features  that 
partially  obscured  lesions,  and  cine- 
graphic  display  further  delineated 
lesions  in  terms  of  time  and  intensity 
of  contrast  change.  As  the  cinegraphic 
format  provided  the  greatest  conspi¬ 
cuity  of  enhancing  tissues,  it  was  used 
to  interactively  guide  ROI  definition 
on  the  time-series  images.  One  super¬ 
ficial  palpable  mass  ( ~  1  cm  in  largest 
diameter)  displayed  signal  intensity 
enhancement,  although  the  enhance¬ 
ment  was  not  distinct  relative  to  sur¬ 
rounding  tissues.  The  mass  region 
was  identified  as  a  capsule  on  the  skin 
that  was  visible  at  MR  ima^ng.  This 
mass  was  occult  at  conventional  x-ray 
mammography  and  US  and  was 
shown  to  be  benign  at  pathologic 
analysis. 

Figure  1  illustrates  data  obtained  in 
a  56-year-old  patient  with  pathologi¬ 
cally  proved  invasive  lobular  carci¬ 
noma  in  two  adjacent  sites  separated 
by  approximately  1  cm.  Mammo¬ 
grams  depicted  only  the  larger  mass. 
An  early-time-point  (ie,  before  con¬ 
trast  intensity  change)  keyhole  image 
(Fig  la)  does  not  depict  the  two  ma¬ 
lignant  foci,  which  were  most  appar¬ 
ent  on  the  late-time-point  subtraction 
keyhole  image  (Fig  lb).  Keyhole  ma¬ 
trix  size  applied  in  this  case  was  256  x 
32  X  16.  Several  keyhole  subtraction 
images  from  the  32-image  time  series 
for  this  section  are  shown  in  Figure 
Ic.  While  the  relatively  rapid  en¬ 
hancement  period  is  visually  appar¬ 
ent  in  the  time  series,  temporal  con¬ 


trast  patterns  of  various  tissues  are 
best  viewed  graphically,  as  shown  in 
Figure  Id.  Mono-exponential  param¬ 
eters  and  fitted  curves  are  also  shown 
in  this  plot.  The  time  constants  of  the 
lesions  (t  =  52  or  55  seconds)  indicate 
that  the  majority  of  signal  intensity 
change  occurs  in  less  than  1  minute. 

The  number  or  location  of  these 
masses  was  not  readily  apparent  on 
precontrast  images  (Fig  la),  although 
their  region  roughly  correlated  with  a 
suspect  mass  seen  on  the  mammo- 
gram. 

An  example  of  a  slowly  enhancing 
lesion  is  Figure  2,  which  was  obtained 
in  a  56-year-old  patient  with  a  prob¬ 
able  fibroadenoma.  Findings  at  MR 
imaging  and  mammography  were 
easily  correlated  for  this  particular 
mass,  which  had  been  sampled  cyto- 
logically  7  years  earlier.  No  change  in 
size  was  seen,  as  noted  in  Materials 
and  Methods.  A  subset  of  the  subtrac¬ 
tion-image  time  series  illustrates 
gradual  enhancement  of  the  lesion 
(Fig  2c).  The  keyhole  matrix  size  ap¬ 
plied  in  this  case  was  also  256  x  32  x 
16.  As  with  the  data  in  Figure  Id,  the 
time-course  curve  is  well  fit  by  an  ex¬ 
ponential  (Fig  2d),  although  for  this 
lesion  the  time  constant  is  substan¬ 
tially  longer  (t  =  119  seconds). 

A  summary  of  fitted  mono-expo¬ 
nential  time  constants  for  a  variety  of 
tissue  types  and  all  patients  is  shown 
in  the  scatter  plot  of  Figure  3.  ROIs 
from  tissues  labeled  ^^normaT'  were 
derived  from  areas  remote  to  the  mass 
that  displayed  sufficient  enhancement 
for  curve  fitting  but  were  otherwise 
normal  in  appearance.  These  were 
typically  regions  of  diffuse  enhance¬ 
ment  from  either  the  ipsilateral  or 
contralateral  breast.  No  more  than 
one  normal-tissue  ROI  was  derived  in 
each  patient.  Two  of  the  five  invasive 
ductal  carcinoma  tumors  were  found 
in  one  patient  (within  one  breast  and 
separated  by  approximately  5  cm), 
and  two  of  the  three  invasive  lobular 
carcinoma  tumors  were  found  in  one 
patient  (Fig  1).  ROIs  of  tissues  adja¬ 
cent  to  invasive  ductal  carcinoma  tu¬ 
mors  were  defined  to  avoid  signal 
blur  from  the  lesion  focus.  As  men¬ 
tioned  in  Materials  and  Methods,  the 
keyhole  matrix  size  was  256  x  16  x  16 
in  three  of  12  patients.  Masses  charac¬ 
terized  in  these  patients  were  three 
invasive  ductal  carcinoma  tumors 
(with  largest  diameters  of  18, 13,  and 
22  mm  and  t  of  14,  20,  and  18  seconds, 
respectively)  and  one  neuroendocrine 
carcinoma  (largest  diameter,  14  mm; 

T,  37  seconds). 

There  was  a  clear  tendency  for  the 
malignant  lesions  to  display  rapid 


contrast  changes  (t  <  60  seconds) 
relative  to  benign  or  normal  tissues. 
However,  one  invasive  lobular  cancer 
displayed  an  exceptionally  long  time 
constant  with  t  of  260  seconds.  Patho¬ 
logic  examination  revealed  the  tumor 
had  a  largest  diameter  of  at  least  1  cm 
(margin  was  positive  at  initial  biopsy), 
gross  size  of  1.2  x  1.0  x  0.8  cm,  no 
angiolymphatic  invasion,  and  no  ex¬ 
tensive  ductal  carcinoma  in  situ;  in 
addition,  17  of  17  lymph  nodes  were 
negative  for  malignancy. 

DISCUSSION 

We  tested  a  technique  with  which 
to  perform  dynamic,  contrast  materi¬ 
al-enhanced  MR  imaging  of  the  en¬ 
tire  volume  of  both  breasts  with  a 
temporal  resolution  of  approximately 
12  seconds.  The  usefulness  of  the 
method  was  explored  in  our  study  of 
12  patients  with  a  total  of  13  breast 
masses.  As  found  in  other  studies,  we 
noted  a  trend  for  malignant  cancers  to 
display  rapid  enhancement  relative  to 
benign  and  normal  tissues,  although 
the  presence  of  one  slowly  enhancing 
cancer  was  observed.  This  would  be 
considered  a  false-negative  case  if  the 
enhancement  time  constant  were 
used  as  the  sole  diagnostic  indicator 
of  malignancy.  The  small  number  of 
patients  in  this  study,  however,  pre¬ 
vents  us  from  drawing  statistically 
meaningful  inferences  regarding  the 
diagnostic  value  of  this  technique. 
Technical  issues  aside,  other  investi¬ 
gators  have  noted  both  false-negative 
and  false-positive  findings  relative  to 
enhancement  rate  (11,12,17,18). 

The  viability  of  this  approach  for 
management  of  breast  disease  neces¬ 
sitates  determination  of  sensitivity 
and  specificity  rates  for  each  applica¬ 
tion.  The  role  that  we  are  currently 
evaluating  is  the  contribution  of  in¬ 
cremental  information  beyond  that 
provided  with  conventional  imaging 
modalities  when  a  suspect  mass  is 
detected.  In  this  capacity,  emphasis  is 
placed  on  characterization  of  previ¬ 
ously  detected  masses  that  are  usually 
more  than  5  mm  in  largest  diameter. 
Lesion  localization  for  eventual  dy¬ 
namic  gadolinium-enhanced  MR  im¬ 
aging  is  hindered  by  nonspecific  le¬ 
sion  properties  in  precontrast  MR 
images  and/or  by  poor  correlation  of 
findings  with  those  of  other  imaging 
modalities,  owing  to  variable  breast 
geometries  and  views.  This  problem  is 
essentially  solved  with  the  3D  key¬ 
hole  approach,  which  has  the  inher¬ 
ent  features  of  complete  coverage 
through  both  breasts  and  high  tempo¬ 
ral  sampling  of  dynamic  contrast 
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Figure  1.  Data  from  a  coronal  image  obtained  through  the  breast  of  a  ^ 
56-year-old  patient  with  two  foci  of  invasive  lobular  carcinoma,  (a)  Early-  ^ 
time-point  (ie,  precontrast)  keyhole  image  does  not  delineate  the  cancer 
tumors,  (b)  Subtraction  image  from  a  late  time  point  (t  =  355  seconds)  in 
dynamic  imaging  shows  the  two  high-contrast  cancer  tumors  {ROWl, 
ROI#2).  (c)  Subset  of  time  series  for  this  section  illustrates  the  relatively 
rapid  contrast  change  in  the  cancer  lesions.  "  =  seconds,  (d)  Time-course 
analysis  of  regions  ROWl  and  ROI#2  with  use  of  mono-exponential 
saturation  model  whereby  both  cancers  are  characterized  by  a  relatively 
short  time  constant  (t  =  52, 55  seconds)  (ie,  a  rapid  contrast  change). 
Images  were  electronically  cropped  to  show  only  the  area  of  interest. 


changes.  These  advantages  would 
also  be  a  benefit  in  other  potential 
applications,  including  use  as  an  alter¬ 
native  breast  screening  modality  in 
problematic  cases  (eg,  dense  breasts), 
as  a  survey  modality  to  search  for  ad¬ 
ditional  malignancies  prior  to  lumpec¬ 
tomy,  as  an  anatomic  guide  to  direct 
biopsies,  and  potentially  as  an  imag¬ 
ing  modality  in  patients  with  sero¬ 
logic  or  genetic  positive  markers  and 
negative  findings  at  mammography. 
Our  encouraging  preliminary  results 
suggest  that  further  investigation 
with  this  technique  is  warranted. 

The  raw-data  domain  of  MR  imag¬ 
ing,  referred  to  as  k  space,  defines  the 
spatial  frequency  distribution  of  an 
object  and  is  a  convenient  format 
with  which  to  describe  various  data 
acquisition  schemes.  In  a  simplistic 
sense,  gross  contrast  features  are  pri¬ 
marily  determined  by  data  near  the 
so-called  center  of  k  space  (ie,  low 
spatial  frequencies),  whereas  features 
with  high  spatial  detail  are  defined  by 
data  from  more  peripheral  k  space. 
These  concepts  are  used  with  great 
success  in  techniques  that  generate 
desirable  contrast  properties  more 
efficiently  than  routine  MR  imaging 
does  (24,25). 

Since  most  changes  in  object  and 
temporal  contrast  are  centralized  in 
k  space,  other  investigators  have  in¬ 
creased  temporal  resolution  in  dy¬ 
namic  gadolinium-enhanced  two- 
dimensional  keyhole  MR  imaging  by 
repeatedly  sampling  only  the  central 


phase-encoded,  ky,  lines  (20).  Images 
with  high  anatomic  quality  and  tem¬ 
poral  resolution  are  produced  by 
combining  dynamic  data  with  previ¬ 
ously  acquired  peripheral  phase-en¬ 
coded  data.  In  application  to  the 
breast,  a  3D  acquisition  is  highly  de¬ 
sirable  for  its  advantages  in  tissue 
coverage,  spatial  resolution,  and  sig- 
nal-to-noise  ratio.  The  second  phase- 
encode  dimension,  kz,  has  an  image¬ 
time  penalty  proportional  to  the 
number  of  sections.  Fortunately,  the 
keyhole  concept  may  also  be  ex¬ 
tended  to  this  dimension  so  rapid  dy¬ 
namic  volume  imaging  is  possible  (21, 
22).  One  version  of  this  sequence  dy¬ 
namically  acquired  only  one-sixteenth 
of  k  space  and  was  applied  in  three  of 
12  patients  in  this  series;  for  the  re¬ 
maining  studies,  one-fourth  of  ky  and 
one-half  of  k^  (ie,  one-eighth  of  k 
space)  were  dynamically  acquired  in 
one-eighth  the  image  time. 


The  main  drawback  to  sampling 
only  a  portion  of  k  space  is  reduced 
spatial  resolution.  This  is  mitigated 
because  dynamically  acquired  key¬ 
hole  data  can  be  complemented  with 
previously  acquired  hill-matrix  data 
to  generate  high-resolution  anatomic 
images.  We  emphasize,  however,  that 
temporally  changing  features  remain 
blurred  by  an  amount  dependent  on 
keyhole  matrix  size  and  lesion  spatial 
dimensions.  That  is,  spatial  resolution 
in  contrast-enhancing  structures  will 
be  limited  by  the  keyhole  resolution 
and  not  by  the  full-matrix  resolution. 
The  rationale  for  splicing  keyhole 
data  into  full-matrix  reference  data  is 
that  it  provides  the  opportunity  to 
view  enhancing  structures  within  a 
high-spatial-resolution  format.  In  this 
format,  morphologic  cues  such  as  stel¬ 
late  patterns  or  irregular  borders  of 
tissues  or  lesions  are  correctly  associ¬ 
ated  with  their  dynamic  enhance- 
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ment  properties.  An  alternative  to  the 
keyhole  method  is  rapid  acquisition 
of  low-spatial-resolution  data  sets  that 
are  reconstructed  and  viewed  inde¬ 
pendently  from  high-resolution  data 
sets.  However,  we  believe  the  com¬ 
plexities  of  keyhole  reconstruction  are 
justified  to  ensure  the  temporal  en¬ 
hancement  profiles  and  spatial  prop¬ 
erties  of  a  region  are  correctly  related. 
Spatial  blur  of  enhancement  that  re¬ 
sults  from  acquisition  of  limited  k 
space  will  spread  along  phase-encode 
directions  to  bordering  tissue.  That  is, 
truly  unenhancing  tissues  can  dem¬ 
onstrate  apparent  enhancement  that 
is  actually  caused  by  bleeding  of  sig¬ 
nal  intensity  from  adjacent  enhancing 
tissues.  Vertical  smearing  of  signal 
intensity  from  small  enhancing  ob¬ 
jects  is  illustrated  in  the  subtraction 
images  of  Figures  lb  and  2b. 

In  addition,  the  apparent  amplitude 
of  enhancement  within  truly  en¬ 
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^  Figure  2.  Data  from  a  coronal  image  obtained  through  the  breast  of 
^56-year-old  patient  with  a  stable  fibroadenoma,  (a)  Early-time-point  (ie, 
precontrast)  keyhole  image,  (b)  Subtraction  image  from  a  late  time  point 
(t  =  355  seconds)  of  dynamic  imaging  shows  contrast  enhancement  in 
the  mass  (ROI).  (c)  Subset  of  time  series  for  this  section  illustrates  a  grad¬ 
ual  contrast  change.  "  =  seconds,  (d)  Time-course  analysis  of  the  mass 
indicates  that  this  fibroadenoma  has  a  moderate  to  long  time  constant 
(t  =  119  seconds). 


The  apparent  enhancement  time 
course  of  the  lesion  is  derived  by  mul¬ 
tiplying  =^(y,z,t)  by  an  ROI  mask  and 
then  performing  spatial  integration. 
The  convolution  theorem  allows  us  to 
relate  this  in  k  space  as 

AROI(£)  =  ROI(0  -  ROI(O) 

=  C(f)/R(ky,k,) 

•  L(ky,k,) 

•  K(ky,k,)dkydk„  (4) 

where  R(ky,kz)  is  the  Fourier  transform 
of  the  ROI  mask.  If  we  asstune  the  spa¬ 
tial  pattern  of  the  ROI  closely  matches 
that  of  the  enhancing  lesion,  then 
R(ky,kz)  «  L(ky  k^);  thus,  the  integrand 
becomes  |  L(y,kz)  p  within  the  keyhole 
area  and  zero  elsewhere.  Finally,  the  full 
matrix  acquisition  is  used  to  normalize 
the  contrast  suppression  effect: 

AROI(0  =  P-C(0-ROI(0); 

where 


hanced  tissue  will  be  diminished.  This 
effect  in  the  phase-encode  dimen¬ 
sions,  y  and  z,  is  demonstrated  with 
the  function  that  characterizes  a  uni¬ 
formly  enhancing  lesion  as 

L(y,z,t)  =  L(y,z,0)  •  [1  +  C{t)] 

•  [1  +  C(f)]},  (2) 

where  C(t)  defines  the  temporal  prop¬ 
erties  with  C(0)  of  0  and  Li(ky,kz,0)  de¬ 
scribes  spatial  properties  of  lesions  in 
k  space.  Changes  in  contrast  in  the 
reconstructed  keyhole  image  of  this 
lesion,  ..^(y,z,t),  are  limited  by  the 
spatial  frequencies  that  lie  within  the 
sampled  keyhole  matrix.  Thus  for 
K(ky  kz)  =  1  within  the  keyhole  range 
and  zero  otherwise,  we  have 

^(y,z,t)  = 

.[1  +  K(ky,k,)-C(0]}.  (3) 


lL(ky,kz)|^dkydk2 


JXeyhole  area 


l^JpuUKspace  1  I ^  ^kydkz  .  (5) 

The  factor  P  represents  the  degree 
of  suppression  in  contrast  change  that 
results  from  sampling  only  a  portion 
of  k  space.  For  the  above  derivation  to 
apply,  the  time  to  complete  each  key¬ 
hole  pass  must  be  short  or  must  be  at 
least  comparable  to  the  period  of 
strong  contrast  change.  The  lesion 
contrast  state  is  considered  nearly 
constant  for  each  keyhole  pass  but  is 
allowed  to  vary  from  pass  to  pass. 
While  coverage  of  full  k  space  may 
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result  in  undersampling  of  temporal 
changes,  such  coverage  is  still  used 
for  normalization  to  provide  a  worst- 
case  bound  for  p.  Within  the  validity 
of  these  approximations,  sampling  of 
partial  k  space  reduces  the  amplitude 
of  the  enhancement  curve,  although 
the  temporal  shape  of  the  curve  is 
correct.  As  a  result,  time-constant  fit¬ 
ting  of  the  shape  of  the  enhancement 
curve  is  relatively  immune  to  the  ef¬ 
fects  of  partial  sampling  of  k  space  as 
long  as  the  sampling  period  is  ad¬ 
equately  short  and  the  signal-to-noise 
ratio  is  sufficiently  high  to  enable  de¬ 
tection  of  attenuated  changes  in  signal. 

For  illustration,  consider  our  de¬ 
fault  settings,  in  which  one-eighth  of 
k  space  is  acquired  with  the  32  ky  x 
16  kz  keyhole  matrix.  3D  acquisition 
with  a  rectangular  30  x  15-cm  field  of 
view  and  4-mm  sections  gives  phase- 
encode  dimensions  of  1.17  x  4.00-mm 
pixels.  The  change  in  suppression  of 
signal  intensity  for  a  single-pixel  (1.17  x 
4.00  mm)  lesion  would  be  p  =  0.13. 
That  is,  if  a  single-pixel  lesion  truly 
changed  contrast  by  100%  (ie,  the  pre¬ 
contrast  level  doubled),  results  with 
this  method  would  indicate  only  a 
13%  change  in  signal  intensity.  For 
larger  lesions,  the  attenuation  effect  is 
lessened  as  the  function  |  L(ky,kz)  p 
becomes  peaked  within  the  keyhole 
area. 

Examples  of  predictions  with  Equa¬ 
tion  (5)  in  larger  lesions  obtained  with 
the  same  acquisition  parameters  are 
p  =  0.46  for  a  4.70  x  4.00-mm  lesion, 
p  =  0.81  for  a  8.20  x  8.0-mm  lesion, 
and  p  =  0.92  for  a  11.70  x  12.00-mm 
lesion.  Composite  cinegraphic  display 
of  anatomic  and  subtracted  images 
partially  overcomes  the  impact  of 
contrast  suppression.  Subtraction  re¬ 
moves  time-invariant  image  features 
that  partially  obscure  lesions,  and 
cinegraphic  display  aids  the  observer 
in  distin^shing  true  contrast  changes 
from  noise  and  other  artifact. 

Data  obtained  in  large  lesions  in 
which  p  ~  1  indicate  signal  intensity 
can  be  enhanced  several  times  more 
than  that  depicted  on  nonenhanced 
images,  but  lesions  with  less  change 
in  signal  intensity  are  also  detectable 
and  well  quantified  with  enhance¬ 
ment-time  constants.  Assuming  true 
changes  in  contrast  are  double  the 
precontrast  levels,  we  anticipate  le¬ 
sions  more  than  5  mm  in  diameter 
(where  p  0.5)  would  be  detectable 
with  this  method  since  the  apparent 
enhancement  would  be  ==50%  or 
more.  Of  course,  weaker  enhance¬ 
ment  or  increased  background  noise 
or  artifact  would  decrease  conspicuity 
of  these  smaller  lesions.' 


The  high  sensitivity  of  contrast- 
enhanced  MR  imaging  for  detection 
of  breast  disease  is  well  known.  Un¬ 
fortunately,  the  modality  has  moder¬ 
ate  specificity  and  is  susceptible  to 
false-positive  findings,  since  some 
benign  tumors  and  tissues  also  exhibit 
enhancement.  Moreover  if  greater 
attention  is  to  be  paid  to  rapidly  en¬ 
hancing  regions,  normal  structures 
need  to  be  appropriately  identified 
that  naturally  exhibit  rapid  signal 
change.  Cinegraphic  images  viewed 
as  individual  sections  (when  acquired 
as  such)  or  drawn  from  a  3D  volume 
set  display  numerous  small  rapidly 
enhancing  regions  representing  blood 
vessels.  To  help  identification  of  ves¬ 
sels  that  can  be  potentially  mistaken 
for  small  lesions,  we  further  exploit 
the  3D  nature  of  the  data  and  pro¬ 
duce  maximum-in  tensity-proj  ection 
images  of  subtraction  images  obtained 
through  the  full  imaged  volume  along 
orthogonal  directions.  Section-to-sec- 
tion  continuity  of  vessels  is  more  eas¬ 
ily  appreciated  in  this  display.  These 
subtraction  maximum-intensity-pro¬ 
jection  images  are  viewed  cinegraphi- 
cally,  thus  they  also  provide  a  quick 
survey  of  the  full  study.  In  this  format, 
changes  in  spatial  and  temporal  con¬ 
trast  of  a  mass  relative  to  other  breast 
tissues,  vessels,  and  blood  in  the  heart 
(usually  partially  included  in  the  im¬ 
aged  volume)  are  readily  apparent. 
Excerpts  of  a  time  series  of  superior- 
to-inferior  maximum-intensity-projec¬ 
tion  images  of  subtraction  images  are 
shown  in  Figure  4a  with  enhance¬ 
ment-time-course  curves  for  regions 
shown  graphically  in  Figure  4b. 

There  are  known  fundamental  tis¬ 
sue  changes  that  provide  the  ratio¬ 
nale  for  the  hypothesis  that  malignant 
breast  cancers  will  preferentially  dis¬ 
play  rapid  changes  in  contrast  over 
most  benign  and  normal  tissues.  Con¬ 
siderable  experimental  data  show  that 
tumors  require  an  enhanced  blood 
supply  to  support  their  rapid  growth 
(26-28).  Thermal  diffusion-based  de¬ 
livery  of  nutrients  sustains  tumor 
nodules  up  to  1-2  mm  in  largest  di¬ 
ameter,  but  their  continued  linear  and 
rapid  exponential  growth  to  larger 
diameters  is  facilitated  by  neovascu¬ 
larity.  Experimental  evidence  that 
tumors  are  angiogenesis  dependent  is 
summarized  in  an  article  by  Folkman 
(29),  in  which  several  studies  are  cited 
that  demonstrate  tumor  size  is  linear, 
slow,  or  stagnant  before  vasculariza¬ 
tion  but  switches  to  rapid  growth  af¬ 
ter  vascularization.  Angiogenesis  does 
not  necessarily  imply  malignancy, 
however,  since  some  benign  tumors 
are  also  angiogenic.  Significant  corre- 
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Figure  3.  Scatter  plot  of  time  constants  de¬ 
rived  from  time-course  data  for  mono-expo¬ 
nential  fits  of  enhancement.  ROIs  were  de¬ 
fined  on  various  tissue  across  12  patients. 
The  N  values  in  the  tissue  type  axis  refer  to 
the  number  of  ROIs  associated  with  each  tis¬ 
sue  type;  only  one  ROI  was  defined  for  each 
discrete  mass.  One  patient  had  two  invasive 
ductal  carcinoma  tumors,  and  another  pa¬ 
tient  had  two  invasive  lobular  carcinoma  tu¬ 
mors.  1  =  invasive  ductal  carcinoma  {n  =  5), 
2  =  metastatic  neuroendocrine  carcinoma 
(n  =  1),  3  =  invasive  lobular  carcinoma  {n  = 
3),  4  =  enhancing  tissue  adjacent  to  invasive 
ductal  carcinoma  (n  =  5),  5  =  fibroadenoma 
(n~l),  6  =  benign  mass  {n  =  3),  7  =  areolar 
or  nipple  tissue  {n  =  8),  5  =  enhancing  nor¬ 
mal  tissue  (n  =  7).  No  more  than  one  normal 
tissue  ROI  was  derived  from  any  given  patient. 


lations  between  tumor  angiogenesis 
and  invasiveness  and  metastatic  po¬ 
tential  of  breast  cancer  have  been 
demonstrated  with  direct  optical  mi¬ 
crovessel  counts  in  breast  tumor 
specimens  (26-28).  There  is  also  some 
evidence  that  precancerous  breast 
tissues  have  a  greater  angiogenic  pro¬ 
pensity,  suggesting  that  angiogenesis 
has  potential  as  a  precancer  marker 
(30).  Factors  that  trigger  and  moder¬ 
ate  angiogenesis  are  not  fully  known, 
but  their  determination  and  potential 
exploitation  for  diagnosis  and  treat¬ 
ment  are  clearly  critical  areas  for  in¬ 
vestigation. 

Evaluation  of  the  clinical  value  of 
imaging  modalities  sensitive  to  breast 
vascularity,  such  as  MR  imaging  and 
Doppler  US  (31),  is  also  an  active  area 
of  investigation.  The  value  of  tempo¬ 
ral  parameters  derived  from  dynamic 
contrast-enhanced  MR  imaging  in 
helping  distinction  of  benign  and  ma¬ 
lignant  masses  is  an  area  of  debate. 
Currently,  there  is  disagreement 
among  investigators  who  have  per¬ 
formed  clinical  studies  to  evaluate  the 
optimal  method  for  performing  con¬ 
trast-enhanced  MR  imaging  for  the 
characterization  of  breast  masses.  As 
mentioned,  technical  factors  related 
to  volume  of  tissue  coverage,  tempo¬ 
ral  sampling  rate,  and  quantification 
of  temporal  parameters  are  variable 
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Figure  4  (a)  Maximum  intensity  projection  of  subtraction  images  in  the  superior-inferior  projection.  Six  nonconbguous  time  femes  are 

shown  from  images  normally  vi^ed  cinegraphically.  ROI#3  delineates  an  invasive  ductal  carcinoma.  Enhancement  m  linear  stoctures  delin 
ea°^vesseTs.“  =  anterr!  1  =  left  P  =  postefer,  R  =  right,  t  =  time  (seconds).  "  =  seconds,  (b)  Corresponding  enhancement-time-course 
for  ROT*;  1-4  illustrate  chanees  in  relative  rate  and  intensity  of  enhancement.  Curves  were  derived  from  ROIs  drawn  on  related  sing 
"cTn  dnSaph";^  (broken  line),  ROm  =  liver  (O),  ROm  =  invasive  ductal  carcinoma  (solid  line), 

ROW4  =  area  adjacent  to  mass  (dashed  line). 


across  these  studies.  The  immediate 
technical  objectives  of  this  study  were 
to  overcome  these  critical  limitations. 
These  objectives  are  largely  met  with 
the  3D  keyhole  technique,  which  of¬ 
fers  full  tissue  coverage  at  high  tem¬ 
poral  sampling  for  subsequent  quantita¬ 
tive  enhancement  time  characterization. 
As  we  have  discovered,  the  ability  to 
dynamically  study  entire  breasts  at  MR 
imaging,  thus  obviating  prospective 
identification  of  one  or  aU  lesions  amidst 
complex  contrasts,  is  a  great  practical 
advantage.  A  drawback  of  this  approach 
is  the  large  data  volume  generated  per 
case.  Currently,  we  resolve  the  data 
management  issue  by  performing  recon¬ 
struction,  display,  and  analysis  on  an 
off-line  workstation  with  customized 
software  tools.  Our  software  has  been 
specifically  designed  to  present  data  in  a 
variety  of  formats  for  specific  purposes, 
such  as  the  survey  of  all  tissues  in  cine- 
graphic  maximum-intensity-projection 
images  of  subtraction  images  or  the  cine- 
graphic  viewing  of  a  single  section  for 
quantitative  analysis  of  individual  tis¬ 
sues.  This  method  also  benefits  from  the 
other  well-known  features  of  breast  MR 
imaging,  such  as  the  ability  to  image 
breasts  with  dense  parenchyma  or  im¬ 
plants  and  to  provide  high-quality  refor¬ 
matted  tomographic  views  that  can  po¬ 
tentially  aid  planning  for  surgery.  ■ 
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